Liquid-liquid phase transition of water in hydrophobic and hydrophilic pores 
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Effect of confinement on the liquid-fiquid transition of water are studied by simulations in the 
Gibbs ensemble. Upon cooling along the liquid-vapor coexistence curve, confined water undergoes 
transition from normal to strongly tetrahedral water via a first order phase transition (as in the bulk) 
or in a continuous way in dependence on pore hydrophilicity. In all cases, transition temperature is 
only slightly shifted by the confinement. This agrees with the experimentally observed weak effect 
of confinement on the temperature of the fragile-to-strong transition of water. 

PACS numbers: 61.20.Ja, 64.70.Ja, 68.08.-p 



The fragile-to-strong dynamic transition (FST) of bulk 
liquid water at about 228 K 1] originates from the change 
of water structure due to the thermodynamic liquid- 
liquid phase transition from normal water to strongly 
tetrahedral water ^ . Crystallization prevents direct ex- 
perimental studies of this phase transition and location 
of its critical point remains unknown. In water mod- 
els with a liquid-liquid critical point at positive pres- 
sure, a sharp change of water properties upon cooling 
at zero pressure corresponds to the crossing the line of 
the heat capacity maxima C^'"^ I^Q- In other models, 
the liquid-liquid critical point is located at negative pres- 
sures and water undergoes a first order phase transition 
upon cooling at zero pressure Water crystalliza- 

tion can be suppressed by confinement and studies of the 
supercooled confined water may clarify the location of 
the liquid-liquid transition of bulk water. Recent exper- 
iments 0, la Ea, H, JJJ show that FST of water in 
various confinements occurs in relatively narrow temper- 
ature interval close to 228 K, where FST is expected for 
bulk water. In the incompletely filled Vycor pores, con- 
fined water undergoes first order liquid-liquid transition 
at 240 K, accompanied by the strong changes of water 
dynamics 7J. In cylindrical silica mesopores, FST occurs 
at about 225 K at zero pressure and its shift with pressure 
follows the behavior expected for the liquid-liquid tran- 
sition i&iE, dOj. In more hydrophobic carbon nanotubes, 
water shows FST at slightly lower temperatures (about 
218 K) HI. Interestingly, the dynamic transition of var- 
ious hydrated biomolecules, accompanied by the onset of 
their function upon heating, occurs approximately at the 
same temperatures (200 to 230 K) jlj|. Obviously, dy- 
namic transition of biomolecules is governed by FST of 
hydration water, which occurs at 222 K and 220 K in the 
case of DNA 1^ and lysozyme 0] molecules, respec- 
tively. 

The phase transitions of fluids are strongly affected by 
confinement and the robustness of the FST temperature 
of water with respect to various confinements seems to 
be surprising. Density functional studies show that the 
liquid-liquid transition of water-like fluid in hydrophobic 
pore is weakly affected by confinement only when the 



fiuid is assume to be homogeneous Computer sim- 
ulations of model water in very narrow pores indicate 
possibility of the liquid-liquid transitions of confined wa- 
ter |l6lll7l |. but liquid-liquid transition of confined model 
water was not found yet. To clarify the effect of con- 
finement on the liquid-liquid phase transition of water, 
we have performed simulation studies of ST2 water |0| , 
whose bulk phase diagram is known in details , in slit- 
like pores of 24 A width. Contrary to the bulk case 5, 4, 
we did not use long-range corrections for Lennard- Jones 
(LJ) intermolecular interactions in pore geometry. The 
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FIG. 1: Phase diagram of water in hydrophobic pore ({7o 
= -0.39 kcal/mol): liquid-vapor coexistence (open circles), 
liquid- liquid coexistence (squares). Liquid- vapor coexistence 
curve of the bulk water is shown by solid circles and solid line. 
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effect of these corrections on the bulk phase diagram of 
ST2 water was found negligible at T > 250 K, whereas it 
causes increase of liquid water density on about 1.5 % at 
lower temperatures. Interaction of smooth pore wall with 
water oxygens was represented by (9-3) LJ potential, 
whose well-depth Uq was varied to reproduce hydropho- 
bic {Uq — -0.39 kcal/mol), moderately hydrophilic {Ua 
~ -1.93 kcal/mol) and strongly hydrophilic {Uq = -3.08 
kcal/mol) pore walls. 
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FIG. 2: Phase diagram of water in moderately hydrophilic 
pore {Uo = -1.93 kcal/mol): liquid- vapor coexistence (open 
circles), liquid- liquid coexistence (squares). Liquid- vapor co- 
existence curve of the bulk water is shown by solid circles and 
solid line. 



of molecules in the two simulation boxes used for liquid- 
vapor and liquid-liquid equilibration was about 850 and 
1600, respectively. The number of successful molecular 
transfers in the course of the simulation runs always ex- 
ceeded 10 N and achieved 100 at high temperatures 
for liquid-vapor equilibration, whereas for liquid-liuqid 
equilibration it was about 2 N. 

The liquid-vapor coexistence curves of ST2 water in 
hydrophobic and moderately hydrophilic pores are shown 
in Figs.l and 2, respectively. In both cases the liquid 
density shows step-like density change by 0.05 to 0.06 g 
cm~'^ at T f» 260 - 270 K. This density step indicates a 
triple point, where a saturated vapor coexists with two 
liquid phases of different densities. A hysteresis observed 
in the temperature range of about 10° clearly evidences 
the first order liquid-liquid phase transition. This finding 
was confirmed by the direct equilibration of the two coex- 
isting liquid phases in the Gibbs ensemble at T > 250 K 
(see squares in Figs.l and 2). Similarly to the bulk case 

Q , the liquid-liquid phase transition of confined water 
shifts to lower pressures (densities) with increasing tem- 
perature and should end at the critical point at some neg- 
ative pressure. The temperatures Tt of the liquid-liquid- 
vapor triple point estimated as a middle of the hysteresis 
loop are shown in Table I. Evidently, the temperature of 
the triple point decreases by just a few degrees due to the 
confinement in hydrophobic and moderately hydrophilic 
pores. Note, that in both pores Tt is essentially above 
the glass transition temperature Tg, estimated from the 
temperature dependence of the heat capacity of liquid 
along the liquid- vapor coexistence curve (see Table I). 

Further strengthening of the pore wall hydrophilicity 
causes appearance of a prewetting transition, which is 
a two-dimensional condensation of a thin water film (of 
about two water layers in the case of nanopores 20] ) . As 
a result, in strongly hydrophilic pore liquid- vapor phase 
transition splits into the prewetting transition and the 
liquid-vapor transition in a pore with a wall already cov- 



The liquid-vapor coexistence of water was simulated 
by direct equilibration of the two coexisting phases in 
the Gibbs ensemble The liquid- liquid phase tran- 

sition was detected by the appearance of a triple point, 
where vapor coexists with two liquid phases of differ- 
ent densities. At the triple point, the liquid branch of 
the liquid-vapor coexistence curve shows a pronounced 
step in density |^ 0|. Additionally, coexistence of two 
liquid phases close to the triple point was studied by 
direct equilibration in the Gibbs ensemble. The use of 
efficient techniques for the molecular transfers (more de- 
tails can be found elsewhere (gj, iSflj) allowed extension 
of the simulated liquid-vapor coexistence curves deeply 
into the supercooled region and direct equilibration of 
two liquid phases at T > 250 K. The total number N 



TABLE I: Temperature Tt of the liquid-liquid-vapor triple 
point and glass transition temperature Tg of bulk and con- 
fined ST2 water. The temperatures T^ — Tt - 30° are ex- 
pected for real water, when the temperature shift of the liquid 
density maximum in ST2 water is taken into account. 
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FIG. 3: Phase diagram of water in strongly hydrophilic pore 
{Uo = -3.08 kcal/mol): liquid-vapor coexistence of water in 
pore interior (open circles), prewetting transition (solid trian- 
gles). Dotted lines show approximately the shape of the coex- 
istence curve expected for 2D transitions. Bulk liquid-vapor 
coexistence curves are shown for ST2 water (solid circles and 
solid line) and ST2RF water (squares). 



ered by about 2 water layers (Fig. 3). The liquid density 
maximum is strongly enhanced in comparison with a bulk 
water and liquid density changes gradually without any 
indications on the existence of a liquid-liquid- vapor triple 
point (Fig. 3, lower panel). Such behavior should be at- 
tributed to the shift of the liquid-liquid critical point to 
higher (positive) pressures. Similar shift of the liquid- 
liquid critical point was observed for the simple van- 
der-Waals model of hydrogen-bonded liquid with the 
strengthening of hydrogen bonds and also for ST2 water 
model with orientational ordering enhanced by the appli- 
cation of the long-range corrections for the intermolecular 
Coulombic interaction (ST2RF model , see lower panel 
in Fig. 3). In the studied strongly hydrophihc pore, wc 
may attribute the shift of the liquid-liquid critical point 
to the strong enhancement of the orientational ordering 
of water molecules near hydrophilic surfaces . Liquid 
branch of the liquid-vapor coexistence curves of ST2RF 
bulk water and of ST2 water in strongly hydrophilic pore 
look like superctitical isobar of the liquid-liquid phase 
transition close to its critical point (Fig. 3, lower panel). 
The inflection of this isobar indicates crossing the line 



of the heat capacity maxima C^"^ emanating from the 
liquid- liquid critical point which is located at about 

280 K for bulk ST2RF water [10. For ST2 water in 
strongly hydrophilic pore G^""^ is located at about 275 
K. So, transition between two liquid phases of confined 
water at zero pressure occurs very close to the temper- 
ature of the liquid-liquid-vapor triple point of bulk ST2 
water (Table I). 

Structural analysis of two liquid phases of bulk ST2 
water, which may coexist with vapor, shows that four- 
coordinated tetrahedrally ordered water molecules dom- 
inates in the low-density phase (phase I), whereas nor- 
mal density water (phase II) is enriched with molecules, 
which have tetrahedrally ordered four nearest neighbors 
and up to 6 molecules in the first coordination shell ^2^ . 
Near hydrophobic surface, these two phases show quite 
different arrangement even in the surface layer (Fig. 4). 
Enhanced tetrahedral ordering in the phase I causes split- 
ting of the water surface layer near hydrophobic sur- 
face (Fig. 4, a). Near the moderately hydrophilic wall, 
the surface layer remains strongly localized and its den- 
sity only slightly decreases upon liquid-liquid transition 
(Fig.4,b). In strongly hydrophilic pore, transition to low- 
density water causes decrease of density in the pore cen- 
ter, whereas the surface layer remains unchanged in a 
wide temperature range (Fig.4,c). So, we may expect 
that the liquid-liquid transition of water near heteroge- 
neous biological surfaces, containing both hydrophobic 
and hydrophilic areas, should affect first of all struc- 
tural and dynamical properties of water near hydropho- 
bic parts of biomolecule. 

Confinement affects strongly the liquid-vapor phase 
transition of ST2 water and its critical temperature de- 
creases by about 50° in strongly hydrophobic and mod- 
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FIG. 4: Density profiles of two coexisting liquid water phases 
in hydrophobic (a, T = 265 K, vertical scale is enhanced) 
and moderately hydrophilic (b, T = 270 K) pores. Density 
profiles of liquid water phases in strongly hydrophilic pore (c, 
T = 300 K: solid line, T = 260 K: dotted line). 
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erately hydrophilic pores (Figs.l and 2) and even more 
in strongly hydrophilic pores (Fig. 3). In contrast, the 
temperature of the liquid-liquid-vapor triple point (or 
the temperature of Cy^) is much less sensitive to the 
confinement. This may be directly related to the large 
slope of the liquid-liquid transition line in temperature- 
pressure plane with respect to zero pressure line. As a 
result, the pressure shift of the liquid-liquid coexistence 
in confinement causes only small shift of Tt or the tem- 
perature of C^°^. This shift is negative in hydropho- 
bic and moderately hydrophilic pore and reflects stronger 
disordering effect of a surface on more ordered strongly 
tetrahedral water phase. The temperature shift is oppo- 
site in strongly hydrophilic confinement, which promotes 
more ordered low-temperature water phase due to the 
strong orientational ordering of water molecules near hy- 
drophilic surfaces. The observed trends agree with higher 
temperature of FST of water in hydrophilic Vycor glass 
(~ 240 0) and sihca pores (~ 225 K '3,'?,'T^) relatively 
to the temperature of FST in less hydrophilic carbon nan- 
otubes (~ 218 K 0). 

Our results indicate, that by tuning the pore hy- 
drophilicity the liquid-liquid critical point may be placed 
exactly at the liquid branch of the liquid-vapor coexis- 
tence curve. Experimental studies of supercooled water 
in pores of various hydrophilicity may clarify, whether 
the temperature of real bulk water singularities (about 
228 K) corresponds to the liquid- liquid- vapor triple point 
or to the distant effect of the liquid-liquid critical point 
located at positive pressure. Note, that the liquid- 
liquid phase transition studied in the present paper is 
the lowest-density one among the multiple liquid-liquid 
transitions of water seen in experiment |23| and simula- 
tions IMS 13- I'^ particular, with increasing pressure, 
ST2 water undergoes transition from normal density wa- 
ter to much less orientationally ordered liquid |22, 123 ■ 
Further studies should clarify effect of this high-pressure 
transition on the properties of confined water and its 
probable relation to the pressure-induced denaturation 
of biomolecules |2^ . 
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